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A model in statistical mechanics, characterised by the corresponding Gibbs measure, is a subset 
of the totality of probability distributions on the phase space. The shape of this subset, i.e., the 
geometry, then plays an important role in statistical analysis of the model. It is known that this 
subset has the structure of a manifold equipped with a Riemannian metric, given by the Fisher 
information matrix. Invariant quantities such as thermodynamic curvature have been studied ex- 
tensively in the literature, although a satisfactory physical interpretation of the geometry has not 
hitherto been established. In this article, we investigate the thermodynamic curvature for one and 
two dimensional Ising models and report the existence of a geometric phase transition associated 
with a change in the signature of the curvature. This transition is of a continuous type, and exists 
for finite systems. The effect may be tested in principle by mesoscopic scale experiments. 
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In the area of parametric statistics, it has long been 
known that a useful and illuminating approach is to view 
statistical models as characterised by differentiable man- 
ifolds, equipped with the Fisher-Rao metric. Of par- 
ticular interest in statistical inference is the notion of 
statistical divergence that measures how separated two 
probability distributions are, which is applied to study 
affinities amongst a given set of populations. It was dis- 
covered by Rao that this divergence can be measured by 
the geodesic distance determined from the Fisher infor- 
mation matrix H. Analogous studies have been made 
in statistical mechanics, and numerous models including 
those exhibiting phase transitions were investigated || 
(for a comprehensive list of references, see S). In par- 
ticular, it was found that, for models exhibiting second 
order phase transitions, the curvature of the thermody- 
namic parameter space Ad diverges at the critical point, 
and that the scaling exponent of the curvature in the 
vicinity of the transition point is identical to that of the 
correlation volume j^]. Furthermore, for simple models 
such as the ideal gas, it has also been shown that the 
geodesic curves on A4 correspond to various equations of 
state for the system j| . 

On the whole, most of the literature on the geometry 
of the thermodynamic parameter space has been devoted 
to analysing properties in the thermodynamic limit, ow- 
ing to the interest in critical phenomena. In this paper, 
however, we investigate finite size (volume) effects in the 
thermodynamic curvature for one and two dimensional 
Ising models. In particular, it is shown that, for a finite 
number (N) of spins, the phase diagram is clearly divided 
into two phases, namely, the positive and negative cur- 
vature phases. Furthermore, this phase separation disap- 
pears as N — > oo. Because the size of the system is finite, 
the transition is continuous and is not accompanied by a 
singularity. Possible experimental realisations, based on 
the idea of generalised equations of state, are considered. 



Before turning to the analysis of finite size effects on 
the thermodynamic curvature for the Ising models, we 
first introduce the basic ideas behind the geometric ap- 
proach to statistical mechanics. For this purpose, we 
consider the van der Waals gas model for vapour liquid 
transitions, because this is a nontrivial model that has 
nonetheless been studied quite extensively. 

In statistical mechanics, we begin by considering the 
probability distribution given by the Gibbs measure, 

p{x\6) = exp^-^ff 4 (x)-lnZ(^ , (1) 

where Hi(x) are Hamiltonian functions on the phase 
space, Z(6) is the partition function, and {9 1 } are ther- 
modynamic variables which may include inverse temper- 
ature, pressure, magnetic field, chemical potential, and 
so on. This density can be mapped to an element in a 
Hilbert space Q by the prescription p(x\9) — > ipo(x) — 
y/p(x\8). Then, because of the normalisation condition 
for p(x\9), the square- integrable function tp$(x) repre- 
sents, for each fixed value of {9 1 }, a point on the positive 
'octant' of the unit sphere S in Hilbert space. By continu- 
ously changing the values of {#*}, this point clearly moves 
inside an r-dimensional subspace of 5. This subspace is 
the thermodynamic parameter space, or, for the above 
Gibbs measure, the maximum entropy manifold Ad, and 
the metric on Ad induced by the underlying spherical 
geometry of S is called the Fisher-Rao metric. In par- 
ticular, if we choose the parametrisation as given in (|l|) 
and write di — djd9 % , then this metric takes the sim- 
ple form Gij = didj \uZ, from which the curvature can 
be computed via standard prescriptions in Riemannian 
geometry. 

Now suppose we consider the van der Waals model, 
where in equation (jl]) we take r — 2, H\{x) the energy, 
H 2 (x) the volume V, and (0\ 9 2 ) = (1/fcT, P/kT). Thus, 
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Z(6) is a function of two variables, namely, temperature 
and pressure. In the case where M. is two-dimensional, 
such as in this example and the Ising models to be con- 
sidered later, the scalar curvature assumes a simple form 
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dflnZ <9i<9 2 lnZ d\\rxZ 
dflnZ dfd 2 lnZ d^lnZ 
dfd 2 lnZ dxdllnZ d\\nZ 



(2) 



where G — det(Gy). 

The phase space integral for the partition function is 
known to give rise to the cubic equation of state for the 
van der Waals model (5) , which includes an essentially un- 
physical region in the parameter space where the pressure 
decreases as the volume is reduced. The conventional ar- 
gument is to apply Maxwell's equal area rule to follow 
through the transition. This is sketched in Fig. 1. Note 
that not all of the region inside the Maxwell boundary 
is entirely unphysical, because in some part the pressure 
increases in decreasing volume. The unphysical region, 
on the other hand, is surrounded by the spinodal curve 
along which d9 2 jdv = 0, where v = v(9 1 , 9 2 ) is the mean 
volume determined by the equation of state. This curve 
also contains the transition point where d9 1 /dv = 0. 




Spinodal curve 



FIG. 1. The equation of state for the van der Waals model. 
The scalar curvature on the parameter space diverges along 
the spinodal boundary which envelopes the unphysical region. 
Thus, one can take the viewpoint that the curvature in some 
sense 'prevents' the entrance to the unphysical domain. 

For the van der Waals model, the scalar curvature TZ 
of the parameter space M. has been calculated , and it 
was shown that the curvature diverges at the transition 
point, as well as along the spinodal curve. This intrigu- 
ing picture can be further analysed by consideration of 
the surface M. on the unit sphere S. In the absence of a 
phase transition, M is a single surface. However, below 
the transition temperature, that is, beyond the spinodal 
curve, the maximum entropy surface M. proliferates into 
three surfaces, associated with the three distinct roots of 
the equation of state. The scalar curvature, in particular, 
is singular along the curve where the proliferation occurs. 
Therefore, suppose we begin with a pure thermal state 
\ip(@)) at high temperature, pure in the sense that it is 



a uniquely defined state in the Hilbert space. Then by 
reducing the temperature, this pure state 'evolves' into a 
mixed state J^j c j\ v j)( v j\ where c 2 determines the prob- 
ability that the system volume is Vj . Then, by a suitable 
measurement to determine the volume of the system, this 
mixed state reduces to a pure state \vk)- This is illus- 
trated in Fig. 2. 
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FIG. 2. A schematic representation of the maximum en- 
tropy surface for the van der Waals model on the unit sphere 
S in Hilbert space. At high temperatures, the surface is 
uniquely defined, while at low temperatures, the surface is 
multi-valued, and can be labelled by the roots of the equa- 
tion of state. Without further information, we cannot identify 
which one of the surfaces the equilibrium state would reach 
beyond the spinodal curve. 

In the case of the van der Waals model, further physical 
arguments in fact rule out two of the unphysical roots for 
the equation of state, and we can single out one physical 
state \vk). However, in a generic scenario of a phase tran- 
sition, in the absence of any symmetry breaking field, the 
pure state characterising the equilibrium state turns into 
a mixed state through a geometric singularity, reflecting 
the multiplicity of the ground state. This is how the 
equilibrium theory for characterising symmetry break- 
ing can be viewed in the geometric framework. Although 
the calculations involved here can be quite intricate, con- 
sideration of thermodynamic limit often leads to useful 
simplifications. 

While the theory of critical phenomena is of great in- 
terest in the approach outlined above, owing to vari- 
ous advances in experimental condensed matter physics, 
there is also a growing interest in understanding the ef- 
fect of a finite system size on thermodynamic quantities. 
This motivates us to study, in particular, such effects 
on the thermodynamic curvature. This will be analysed 
here by consideration of the Ising models in one and two 
dimensions. Let us first study the one-dimensional Ising 
chain. The Hamiltonian for the system is given by 
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■ ph. = /3 ai<Tj + i + h ai 



(3) 



where {cr, = ±1} are the spin variables, f3 = l/(fcT), 
and the two-dimensional manifold A4 is parametrised by 
(9 1 , 9 2 ) = (f3, h). The components of the Fisher-Rao met- 
ric per spin are obtained by differentiating lnZ(/3, h)/N: 
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Gij = —didj {In [(cosh/i + + (cosh h - r,) N ] } , (4) 



where r] — y sinh 2 h + e~ 4/3 . Here, we set Boltzmann's 
constant k — 1. If we compute the metric in the ther- 
modynamic limit N — > oo, then the resulting expression 
simplifies, and we obtain the thermodynamic curvature 
H , given by 1Z = 1 + ?/~ 1 cosh h, which is always positive. 
However, if we consider finite size effects, we observe that 
the curvature is no longer strictly positive. In particular, 
it was noted in || that for a given point on parameter 
space, 1Z could decrease and eventually become negative 
as N was reduced. Indeed the transition to large nega- 
tive values over a small range in N is quite marked as 
one may observe in Fig. 3 where we consider the point 
{T = 0.3, h = 0.4}. 



bination of the products of central moments for terms 
in the Hamiltonian, and thus calculation of 7Z is quite 
amenable to Monte Carlo analysis. 




FIG. 3. A plot of thermodynamic curvature 1Z for the 
one- dimensional Ising model, with {T = 0.3, h = 0.4}. If 
the number (N) of spins is reduced from 35, the curvature 
rapidly grows to large negative values. When N is increased, 
1Z approaches its thermodynamic limit 1Z ~ 3.63. 

The result in Fig. 3 suggests the existence of a phase 
separation in terms of positive and negative curvatures. 
Indeed, because the components of the metric can be ob- 
tained explicitly for arbitrary finite N 0], we can com- 
pute the curvature 1Z as a function of N and then nu- 
merically extrapolate the phase boundary along which 
the curvature vanishes. The result is shown in Fig. 4 for 
several values of N . Further numerical analysis shows 
that the parameter space is indeed separated into these 
two phases. Before seeking any physical consequences of 
this geometric phase separation, we would like to anal- 
yse the case of the two-dimensional Ising model, in order 
to indicate the universality of the phenomenon at least 
within the context of ferromagnetic spin models. 

In the two-dimensional case, we also consider the 
Hamiltonian given in (^|) parametrised by the inverse 
temperature and magnetic field, except the spin variables 
{<Ji\ are now defined on a toroidal lattice. For the two- 
dimensional Ising model with an applied field, an ana- 
lytical expression for In Z is not available for sufficiently 
large N. However, a Monte Carlo simulation turns out to 
be quite effective in analysing problems of this kind ut . 
In particular, if we return to the expression given in (^), 
we notice that each entry in the determinant is a com- 
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FIG. 4. The geometric phase diagram for the Ising chain 
in the (T, h) plane. The curvature is positive in the high 
temperature phase, negative in the low temperature phase, 
and vanishes along the curves indicating the phase boundary 
for N = 8, 12, 16. 

We simulated the Ising spins on 8 x 8 and 16 x 16 
lattices, using a standard Metropolis algorithm, and the 
results were combined from a series of runs using alter- 
nately sequential and random site updates, and hot and 
cold starts. The values of 1Z were computed for a range of 
values of (T, h) , and the results clearly indicate the pres- 
ence of a transition between positive and negative cur- 
vature regions. Numerical values for the curvature near 
this phase boundary, including statistical errors from a 
series of 10 runs with the differing initial conditions and 
update procedures mentioned above, are presented in Ta- 
ble I. The relatively large errors are a result of significant 
cancellations between the third order cumulants used to 
construct the curvature. An interesting point to note 
is that the boundary is fairly stable for the two system 
sizes, in contrast to the shift observed in the ID case. 
This is presumably related to the presence of the criti- 
cal point at (T c ~ 2.27, h — 0). Preliminary results for 
a similar analysis of the 3D Ising model also indicate 
a clear presence of a negative curvature region for low 
temperatures, although this region in parameter space is 
somewhat smaller than for the 2D case above. 

In the foregoing analysis we have demonstrated the 
existence of positive and negative curvature phases for 
one and two dimensional Ising models. Whereas the dis- 
continuities which are manifested in conventional phase 
transitions can only occur in the thermodynamic limit, 
these qualitative geometrical transitions occur even in 
the simplest finite-dimensional models. 

The important issue we now address is whether there 
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is any physical or observational consequence of this tran- 
sition. While the answer appears to be yes insofar as 
the curvature is at least indirectly measurable, a more 
intriguing consequence arises once we recall the fact that 
the deviation of a collection of nearby geodesic curves is 
conditioned by the Ricmann curvature tensor. In partic- 
ular, if we consider a set of nearby geodesies in a region 
of A4, then these geodesies will tend to deviate apart if 
sgn(7£) < 0, while the converse is true if sgn(7£) > 0. Fur- 
thermore, the geodesies on M. can be viewed to play 
the role of generalised equations of state corresponding 
to different initial conditions. This can be seen from the 
fact that the leading term in the Taylor series expan- 
sion of the relative entropy is given by the infinitesimal 
line element ds 2 = GijdB l dBK In other words, we have 
S(p(9)\p(9 + d9)) = \ds 2 in the limit dO -> 0. Hence 
we see that a geodesic joining two nearby points on A4 
corresponds to a trajectory that extremises the entropy, 
and thus has a natural parallel with an equation of state. 

On the other hand, for any given point in Ai there are 
infinitely many geodesic curves that pass this point, each 
of which is subject to an initial condition, which may 
include the requirement of an adiabatic or isothermal 
change of the system. In this case, two nearby equations 
of state would deviate away if the sign of TZ is negative. 
Therefore, an experimental verification of the phase sepa- 
ration requires an extensive analysis of equations of state 
for various different processes. Given such data it may 
be possible to identify qualitatively different behaviours 
for the two phases. 



8x8 


h=0.02 


0.18 


0.34 


0.50 


0.66 


T=2.2 


-2.8 


-9.2 


-4.5 


-5.8 


-0.1 


2.3 


-1.9 


-5.9 


-4.7 


-2.5 


+1.5 


2.4 


+0.6 


-2.8 


-1.2 


-1.4 


+0.1 


2.5 


+1.1 


-2.8 


-1.1 


-1.3 


+0.7 


2.6 


+2.8 


-2.4 


-1.3 


-0.5 


+3.3 


2.7 


+4.1 


-2.5 


-0.7 


-0.2 


+0.2 


2.8 


+5.6 


-1.2 


-0.1 


+0.1 


+1.1 


2.9 


+7.0 


-0.7 


+0.7 


+1.0 


+2.0 


3.0 


+7.3 


+2.4 


+0.1 


+ 1.5 


+3.0 



16 x 16 


h=0.02 


0.18 


0.34 


0.50 


0.66 


T=2.2 


-58.3 


-7.8 


-5.7 


-2.8 


-2.5 


2.3 


-55.1 


-7.1 


-5.6 


-3.5 


-2.6 


2.4 


-32.4 


-6.4 


-5.3 


-2.4 


-2.1 


2.5 


-6.9 


-2.8 


-2.6 


-1.2 


+0.2 


2.6 


+13.2 


-2.1 


-0.7 


-0.9 


+1.1 


2.7 


+22.1 


+0.9 


-0.8 


-0.1 


+1.6 


2.8 


+24.3 


+1.8 


-1.0 


-0.1 


+2.5 


2.9 


+25.8 


+2.4 


+0.6 


+0.1 


+2.6 


3.0 


+23.9 


+2.9 


+0.7 


+0.3 


+1.3 



TABLE I. The values of the curvature TZ (±0.6) of the 
parameter space for the two-dimensional 8x8 and 16 x 16 
Ising lattice. The results clearly indicate the phase separation, 
with a negative curvature phase for low temperatures. 



In summary, we have shown the existence of phase sep- 
arations in terms of the signature of the thermodynamic 
curvature for finite Ising spins in one and two dimensions. 
Owing to the fact that TZ is expressible in terms of corre- 
lation functions of order < 6 (< 3 in the present case), a 
Monte Carlo computation is straightforward, in contrast 
to the renormalisation group approach which generally 
requires knowledge of an infinite number of correlation 
functions. Therefore, it would be of interest also to anal- 
yse the curvature for lattice gauge theories in a finite 
volume. Finally, we have argued that geodesies on M 
generalise the notion of an equation of state, for arbitrary 
initial conditions, and thus experimental realisation may 
be obtained in principle by the analysis of equations of 
state for these systems at mesoscopic scales. 

Another interesting open question in this regard is the 
direct physical interpretation of the curvature TZ. We 
note that TZ is determined by the parallel transport of 
tangent vectors over A4. In particular, by following a 
closed contour C on M, a tangent vector is rotated by 
an angle given by the integral of TZ over the area inside 
C (holonomy). In ordinary thermodynamics, Maxwell's 
relations guarantee that any closed contour C is ther- 
modynamically trivial, provided C does not surround a 
critical point or straddle the spinodal curve. Therefore, 
if there is a physical interpretation of tangent vectors 
over the maximum entropy surfaces, then the geometric 
notion of holonomy leads to a new thermodynamic prop- 
erty whereby a system does not return to its original 
state after going around a contour, a concept analogous 
to magnetic hysteresis. 
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